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ABSTRACT 



The spacial distribution of fields and currents in confining theories can give direct 
evidence of dual superconductivity. We would like to discuss the techniques for 
finding these properties and calculating the superconductivity parameters in 
lattice simulations. We have seen dual Abrikosov vortices directly in pure U(l) 
and SU(2) and others have also seen them in SU(3). In the non-Abelian cases 
' the system appears to be near the borderline between type I and II. We also 

discuss the response of the supercurrents to external fields. 
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qq ' 1. Introduction 

Lattice gauge theory offers the prospect of exploring dual superconductivity^ in 
depth as a confining mechanism. The existence of a dual Abrikosov vortex between a 
static quark-antiquark pair leads to a linearly rising potential energy between them 
and hence confinement. These have been seen in U(l), SU(2) and SU(3) lattice gauge 
simulations thereby supporting this picture of confinement. 



In ordinary superconductivity, the primary issue is the spontaneous breaking of 
the electromagnetic U(l) gauge symmetry (SSB) signaled by the non-vanishing of the 
vacuum expectation value of a charged field. An immediate consequence is that the 
curl of the vector potential is proportional to the curl of the electric current known 
as the London relation. The London relation is violated only near the boundaries 
of superconducting material within a distance of the order of the coherence length. 
Combining the London relation with Maxwell's equations gives mass to the electro- 
magnetic field. The Meissner effect and infinite conductivity follow. 

The theory of superconductivity entails (i) the identification of what symmetry is 
broken and what are the relevant coordinates, (ii) the mechanism that leads to the 
the instability and hence SSB, e.g. BCS theory, and then (iii) an effective theory of 
the currents and fields in the broken phase, e.g. Ginzburg-Landau [GL] theory, which 
predicts the spatial consequences of the broken symmetry!. 

Our goals are similar to these. We hope to identify relevant coordinates and a 
corresponding effective theory. Experience has shown that an analytic approach is 
very difficult. The Villain form of the U(l) action yields some analytic results^, but the 
Wilson form for U(l) and SU(N) has not. Lattice simulations are capable of building 
evidence for or against this picture and perhaps even be a guide for analytic efforts. 



The dual London relation found in pure gauge U(l)S and SU(N)li simulations gives 
evidence that an effective lattice dual Higgs theory is a credible candidate mechanism 
for confinement. 

Fields and currents exist only near the boundaries of dual superconductors over a 
distance scale set by the London penetration depth \ d . It is these spatially transient 
effects that provide spatial structures where the local properties of dual supercon- 
ductivity can be studied and that is the focus of this work. A dual Abrikosov vortex 
is an example of such a structure. The core is 'normal material' and the flux tube 
connecting quark and antiquark is a spacially transient boundary effect. 

2. Abrikosov vortices in the lattice Higgs effective theory 

I would like to review briefly the salient features of the lattice Higgs theory in 
the tree approximation. It is relevant to note that in this model, the London relation 
is a direct consequence of the spontaneous breaking of the U(l) gauge theory. We 
compare the simulations of pure gauge theories with the dual of this model. 

Consider the lattice action in a standard notation: 

S = Yl (! - cos W^)) - KY.ifWe^Mz + + H.c.) + £ V Higgs {\(f>{x)\ 2 ). 

X,fJ,>U X,fl X 

(1) 

We are interested how the fields and currents respond to a magnetic monopole anti- 
monopole pair. We obtain a classical solution by minimizing the action. Our initial 
configuration contains a single closed magnetic monopole loop in the x^,X4 plane 
analogous to the Wilson loop projector in the pure gauge simulations. Our algorithm 
then rejects all updated links that change the initial monopole configurations. We use 
the method of simulated annealing, slowly increasing /3, holding A 2 /a 2 = j3/n = l/e 2 K 
constant, where A/a is the London penetration depth in lattice units and a is the 
lattice spacing. 

The electric current is given by 

n 3 

— Jt(x) = ImU*(x)e ie ^ x) (t)(x + e M )). (2) 
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Writing the Higgs field 4>{x) = p(x)e lu( * x \ the curl of the current is hence given by 

^(A,r u (x) - a u j;(x)) = 

p(x)p(x + e (M) ) sm[-w(x) + 6^(x) + u(x + e ( ^)] + 

p(x + e M )p{x + e M + e M ) sm[-u{x + e M ) + 6 v {x + e M ) + u{x + e M + e M )] - 
p(x + e^)p(x + e W + e (*0) sm[-uj{x + e M ) + 6^x + e M ) + u{x + + e (i/) )] - 
p(x)p(x + e M ) sin[-cu(x) + 9 v {x) + uj(x + e M )], (3) 



where is the lattice derivative. Compare this with the electromagnetic field tensor 

ea 2 F^ u = srnlO^x) + 9 v (x + e<">) - 9„{x + e^) - 9 u (x)). (4) 

If the U(l) gauge symmetry is spontaneously broken, these two quantities are equal 
which is the London relation. To be more precise: if (i) p is nonvanishing and inde- 
pendent of position (absorb the normalization into k) and (ii) 

sin[0 + 2rm] « 9, (5) 

then 

J> = - ^(A^(x) - A v J e (x)) = —- 2 = ne m - (6) 

where J r fiu a 2 is called the fluxoid and e m is the Dirac monopole charge. Condition (i) 
is satisfied by the SSB potential 

VWl0(x)| 2 )^|0(x)| 2 = l, (7) 

and (ii) is satisfied if the lattice is a good approximation, yet still allowing singular 
configurations when n / giving the quantized fluxoid. (A Dirac string would act 
differently, adding 2irn to both quantities.) A 'mexican hat potential' can also give the 
London relation but it would be violated if the plaquette is within the coherence length 
of a superconducting-normal boundary where the Higgs field has a non-vanishing 
gradient. 

Eqn|| together with Maxwell's equations, gives the Meissner effect B — X 2 V 2 B = 
0, where A 2 = m~ 2 \ infinite conductivity E = X 2 A^J (assuming p = 0); and an 
Abrikosov vortex 

B z -\ 2 (Axf) z = n e -^5 2 x± , . (8) 

Fig.0 shows the profile of the R.H.S. of eqn.([3]) and eqn.(^j) in the directions 
perpendicular to the 5x5 magnetic monopole loop on a 12 4 lattice, with /3/k — 1. 
We used the constrained form of the Higgs potential, eqn.(^), which corresponds to 
an extreme type II superconductor. The graphs show the expected behavior: i.e. the 
equality of the two quantities everywhere except at r = where they should differ 
by 27r. There are significant violations only at r = and 1 due to the breakdown of 
Eqn. @. 

Lattice artifacts can violate the London and fluxoid quantization relations. But 
ignoring these controllable effects, the interior surface spanned by the monopole loop 
is in the normal phase since the London relation is violated there. All other regions are 
superconducting. This translates in the following sections in which the Wilson loop 
projects out the normal phase in the plane spanned by the loop and the remaining 
regions are a dual superconductor. 

3. Dual superconductivity in pure U(l) gauge theory 



We now turn to the dynamical simulations in pure U(l) gauge theory given by the 



Fig. 1. Profile of ea 2 F\2(x), eqn.(Q), and — (a 4 /eK)(Ai Jf (x) — A2Jf(x)), eqn.(||), as a function of 
distance perpendicular to the plane of the monopole current loop in lattice units. 



first term in the action, eqn.(|l|). We have given evidence in ref.Q that a dual London 
relation is satisfied: 

I ft f 

Q m = -e^ T F aT - %{ApJ?{x) - A„J?(x)) = -^5 X ^ 5 Xvfi . (9) 

Magnetic monopoles are the current- carriers responsible for the dual superconductiv- 
ity. These monopoles are defined in a 3-volume by the DeGrand-ToussaintB construc- 
tion, separating multiples of 2ti from the plaquette angle, (— n < O^x) < n): 

6 ia/ (x) = e iiV {x)+2Tm liV (x)\ J™{x)=e^ T {6 aT {x + e {u) )-Q UT {x)) (10) 

We associate JJ?(x) with a link on the dual lattice, making world lines which define a 
conserved current density. The key idea is to measure the line integral of J™ around 
a dual plaquette, giving a 2 (curlJ m )^ u . The correlation of this with a Wilson loop, 
gives a signal for the solenoidal behavior of the currents surrounding the electric flux 
between oppositely charged particles. 




Fig. 2. Operators for (a) ea 2 F 34 (x), and (b) (a 4 /e m )(Ai J^(x) - A 2 Jf(x)). 

Fig. |2| shows the lattice operators for the electric field and the curl of the magnetic 
monopole current. The longitudinal electric field, (a), is given by a z,t plaquette 
which is depicted by a bold line for fixed time. The curl of the magnetic monopole 



Fig. 3. Profile of ea^F^x) and — (a 4 /e m )(Ai J™(x) — A2 J™(x)) as a function of distance perpen- 
dicular to the plane of the Wilson loop in lattice units. 

current, (b), is built from four 3- volumes which appear as squares since the time 
dimension is not shown. Passing through the center of each square is the link dual 
to the 3-volume. One takes the monopole number n in each 3-volume and associates 
the value ne m with the corresponding dual link. The 'line integral' around the dual 
plaquette completes the picture. Notice from this construction that E and A x J m 
take values at the same location within the unit cell of the lattice, both are indicated 
by the bold face line in the z direction. Both operators are in defined on the same 
z, t plaquette. 

Fig.^ shows the profile for the quantities making up the London relation on the 
same plot as for the effective lattice Higgs case, Fig.[l]. We did 800 measurements on 
a 12 4 lattice, (3 = 0.95, and a 3 x 3 Wilson loop was used to project onto the q, q 
sector using methods described in Ref.(S). We did a \ 2 fit °f the London relation, 
eqn.(||), using the r ^ points and determined = 0.49(3). that determined the 
coefficient of the delta function in eqn.@, ne = 1.07(7), compared to the expected 
value = 1.03. All these features correspond directly to the classical solution of 

the lattice Higgs model reinterpreted as a dual theory. 

4. Generalization to pure SU(2)and SU(3) gauge theories 

In ref.(@) we applied these same techniques and showed that dual Abrikosov vor- 
tices also occur in SU(2) pure gauge theory in the maximal Abelian gauge. Matsubara 
et.al.i confirmed these results with better statistics and generalized them to SU(3). 
The SU(2) link matrices are U^(x) and the action is 



S = P E (1-krEW*)). 



(11) 



The maximal Abelian gauge is defined by maximizing the quantity 



R = Y,Tr[a 3 U,(x)a 3 Ul(x)}. 



(12) 



Fig. 4. Same plot as in fig.|| but for SU(2), in the confined phase 



Fig. 5. Same plot as in fig.[||4]but for SU(2), in the deconfined phase. 

The Abelian link angle is then taken as the phase of [t^ M (^)]n and the calculation can 
proceed with little changed. 

Recently we have done the calculation at finite temperature in order to check this 
picture on each side of the deconfining phase transition! which we report here. The 
results are shown in the confining phase, fig.|4], (3 = 2.28 and the deconfining phase, 
fig|| (3 = 2.40 on a lattice 4 x 17 2 x 19, with 800 measurements for each case. Gauge 
fixing required about 600 sweeps for each configuration. Fig.[| shows an important 
difference from the U(l) case. Whereas in the U(l) case a linear combination of 
E and curlJ m for r / can vanish giving the London relation, it is clearly not 
possible here. The behavior of —curlJ m does not match that of E. We interpret 
this discrepancy as a signal of a non-zero Ginzburg-Landau coherence length, 
Whereas all evidence in the U(l) case points to the extreme type II limit, i.e. the 
superconducting order parameter turning on at the surface, the evidence here is that 
the order parameter turns on with a behavior tanh(x±/^d) The value of the coherence 
length is approximately the radius where the London relation is restored. 

The simulation shown in fig.|5| was identical to fig.|] except that (3 was increased to 
2.40 putting us in the deconfining phase. The dramatic decrease of curlJ m results in 
the failure of the dual superconductor interpretation as expected for the deconfining 
phase. Also E falls more slowly with radius. 

The behavior of all the SU(2) and SU(3) examples is similar to that shown in fig.|4]. 



The interesting conclusion is that = \d/£,d ~ 1 (no relation to Higgs k) For a type 
II dual superconductor n d > l/y/2 and a type I otherwise. These simulations indicate 
that the non-Abelian dual superconductors lie at the borderline between type I and 
type II. However a more systematic study is needed to pin this down. 



5. External Electric Field 



Many interesting superconducting properties can be elucidated by studying the 
properties of the material in the presence of a magnetic field. We report a preliminary 
look at the corresponding problem of dual superconductivity in the presence of an 
external electric field. The Wilson loop provides such a field by projecting a qq out 
of vacuum configurations. But it may be interesting to see the spontaneous breaking 
of translation invariance as the electric flux forms a vortex rather than impose the 
vortex position at a given location and with quantized flux. Type I and type II dual 
superconductors respond very differently to background uniform external fields. 

For a periodic U(l) lattice the sum of the plaquette angles over any plane is 
identically zero. A classical uniform electric field for U(l) on a particular time slice 
can still be obtained by constraining one plaquette angle in each z, t plane to the 
value (l/(N z N t ) — 1)6 C - Then the remaining (N z N t — 1) plaquette angles will take the 
common value 8 C / (N z N t ) = ea^F^ giving a uniform electric field on all but one time 
slice. (N labels the lattice dimensions.) Choosing 9 C = it gives the largest field. For an 
8 4 lattice and (3 = 1 for example, there would be enough electric flux to form about 3 
vortices in the x, y plane in the dual superconducting phase. This field configuration 
can also be obtained by multiplying the plaquette that was singled out in the action 
by a minus sign and such configurations have been studied in non-Abelian theories^. 
We avoided an alternative method of imposing an external field by introducing a 
non-zero equilibrium value locally for each plaquette angle since we eventually want 
to see the field break translation invariance. Turning on interactions for (3 < 1 brings 
up other interesting featured 

Our goal here is to try to see a signal showing that curlJ m responds to the external 
field. The immediate problem is that the sum of (curlJ m ) xy over the any x, y plane 
is identically zero, and the sum of E z is not. Yet we expect the London relation to 
be satisfied. The only possibility is that translation invariance is broken which is 
expected since vortices segregate the superconducting and normal phases. 

We make the following rough ansatz that the local London relation is due to a the 
alignment of local current loops in the external field. This suggests that we truncate 
curlJ m to include only values ±2, ±3, ±4 representing the winding of the current 
around the dual plaquette. We denote this J™ (2 + 3 + 4). 

Now we get a large signal for curl J™ (2 + 3 + 4) as shown in fig.|] with a sign 
that agrees with the sign in fig.|3| for the superconducting region, r ^ 0. Further if 
we recalculate fig.H with the truncated current, we find a large suppression at r = 



Fig. 6. Average of ea 2 F 34 and -(a 4 /e m )(Ai J 2 m (2 + 3 + 4) - A 2 J{ n (2 + 3 + 4)) on each time slice 
for an external field = 7r/64 (the horizontal line). The constrained plaquette is at t = 1, the field is 
classical at t = 2, 8 and (3 — .99 for t = 3 — 7. 

and a moderate enhancement at the other points. In other words this choice biases 
in favor of the dual superconducting phase. 
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